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ABSTRACT. We consider a porous medium flow in which the gas is initially 
distributed in the exterior of an empty region (a hole) and study the final stage 
of the hole-filling process. From general theory it is known that this hole- 
filling is asymptotically described by a self-similar solution which depends on 
a constant determined by the initial configuration. In general, this constant 
must be found either experimentally or numerically. Here we give an example 
of a one-dimensional symmetric flow where the appropriate constant is obtained 
explicitly. 

PACS number(s): 47.56+r 

Self-similarity arises in a natural way in the description of critical behavior 
in various physical problems. Many examples can be found in the book of G. I. 
Barenblatt [1]. In this note we consider the focusing problem for porous medium 
flow. In this problem the material is initially distributed in the exterior of an 
empty region (a hole) and one is interested in the details of the final stage of 
the hole-filling process. For axially symmetric flows it is known [2] that this 
process is ultimately self-similar. Here we present an explicit calculation of the 
development of self-similarity in the focusing of a particular one-dimensional 
flow. 

The evolution of the (scaled) density U(r, t) of an ideal gas flowing in an 
axially symmetric homogeneous porous medium is governed by the degenerate 
nonlinear diffusion equation 

d t u = (dl + — -d r ){u m ), 

r 

where d is the dimension of the space and m is a constant derived from the ideal 
gas law. Using the ideal gas law we can replace the density U by the (scaled) 
pressure 

V = -^—U m ~ l 
m — 1 

whose evolution is governed by the equation 

d t V = (m - l)V(3 2 rr + — 5 r )V + {d r Vf. (1) 

r 

In the context of gas flow m > 2, but other values of m occur in various ap- 
plications. Here we will assume that m > 1 which assures that the speed of 
propagation of disturbances from rest is finite (slow diffusion). For m < 1 the 
corresponding speed is infinite (fast diffusion) [3] . 

In the focusing problem for (1) we seek a solution in [0, oo) x (r, oo) for some 
t e R such that 

F(-,r) = Vb(-)on[0,oo), (2) 
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where Vq is a given function which satisfies 

V(r\i = 0on [0,a]U[6,oo) 
Vo{r >{ >0on( a ,6) 

for some < a < b < oo. It is known [3] that problem (1),(2) possesses a unique 
continuous generalized solution. As time increases from t = r material flows 
outward from r = b and inward from r = a. There is a non-increasing inner 
interface curve r = a(t) and a non-decreasing outer interface curve r = b(t) with 
o(r) = a, 6(r) = 6, and 



V(r,t) 



= on [0,a(t)] U oo) 
> on (a(t),b(t)) 



Both interface curves become monotone in finite time and there is a finite T > r 
such that a(t) > for t < T and a(T) = 0. 

For normalization we will assume that T = so that r < 0. It is shown in 
[4] that there exists a one-parameter family {g c (r, t)} of self-similar solutions 
to (1) defined for c G R + and (r,t) G [0, oo) x (— oo, 0] which focus at t = 
(the Graveleau solutions). Specifically, there exist numbers a*(d,m) G [1,2) 
and j(d, m) G such that for t < 



flcM) = — <p(c??), (3) 



where 
and 

g c satisfies 
where 



v5 



5c (r, £) 



> on (7, 0) 
= on (— oo, 7) 

7] = tr~ a \ 

= for r G [0,p c (t)] 
> for r G (p c (i),oo) 



is the interface. The function <p and the similarity exponent a* are obtained by 
solving a nonlinear eigenvalue problem and, in general, must be found numeri- 
cally. However, in one space dimension (d = 1) we have a*(l,m) = —7(1,771) = 
1 and each Graveleau solution is a pair of converging plane waves given by 

g c (\x\,t)=c{\x\+ct}+, (4) 

where {•}+ = max(0, •). Thus, for d = 1, 

77 = R and ^ (c) = _c(1 + c) - (5) 
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The Graveleau solutions describe the asymptotics of the focusing problem 
(1),(2) in the following sense [2]. 

Let V denote the generalized solution of the initial value problem (1),(2). 
There exists a c* G R + (depending only on d,m 7 and Vo) such that the inner 
interface 

a(t) ~ p c ,(t)as t SO. 

For each fixed i] e (— oo, 0] the profile V approaches the profile g c * as r \ 0, 
i.e., 

V(r,nr a ) ~ g c *(r,nr a )as r \ 0. (6) 

The constant c* depends on the initial function Vo and so must be determined 
on a case-by-case basis, usually numerically. However, in one particular case we 
can determine c* explicitly. Consider the focusing problem for d = 1 where we 
start at some time r < with symmetrically placed point masses at x = ±£ 
and assume that focusing takes place at T — 0. In view of the symmetry it will 
suffice to consider x > 0. 

The solution to (1) corresponding to a mass M initially concentrated at 
(£, t) is given by 

where 

m - 1 2m(m +1) V B 

M = ,4(m+l)/2(m-l) B -l/2 f ( cos 0) (m+l)/(m-l )de _ 

^0 

Note that Vm{x, t, r) > only on the interval 

f-iJ(t-T) < a; < ^ + i2(t-T). (7) 

This solution is due to Barenblatt [5] and is a self-similar solution of the first kind 
meaning essentially that the similarity exponent can be found by dimensional 
analysis. In contrast the Graveleau solution is a self-similar solution of the 
second kind whose similarity exponent a* cannot be obtained by dimensional 
analysis alone [1]. 

In the focusing problem we require that the left hand boundary of the sup- 
port of Vm just reaches the origin at t = 0, i.e., that 

i = R(-r) = ^(-rr. (8) 
For given (£, r) this condition determines the admissable mass since 

A= _^e_ 
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Using (8) we can rewrite Vm as 



2(-T)V r) \ (l-(/ T )W/ + - 

Note that, in view of (7), Vm{x — £, t — r) — > as x — > for any f 6 (r, 0]. 
If we introduce the dimensionless quantities 



{-t)Vm t x 

-^' ni = 7' n2 = ? 



we can rewrite (9) in the form 



n-^.iw.l<i-no^{i-£^} + . 

where II — > as II 2 — > for each III . Set 

n* = ni/n 2 . 

Then 

n = I<i-nn W -{ 1 - TI <I^} + . 

With II* held fixed, it follows from Taylor's theorem that 

n = n 2 (i - /3iT) + oijij) as n 2 -> o. 

Therefore 

$(ni, n 2 ) - n 2 (l - /?IT) as n 2 -> with IT constant. (10) 
To interpret this result in the original dimensional variables we set 

n = — and p = (11) 

X —T 

Note that r\ < 0,p > 0, and II* = -pry. It follows from (10) and (11) that for 
each fixed rj 

V M {x,t;,t,T) ~p/?a;{l +pprj} + as x — > 0. 
Using (5) we can rewrite this as 

V M (x,£,t,T) ~ —ip(p0v) = 9 P p(\x\,t). 
Thus (6) holds with 

c* = p(3. 

Similarly, for the inner interface we have II* <~ 1/(3 so that 

a{t) ~ pfi(-t) = p pfl (t). 



Two aspects of the explicit computation of c* are essential. Since we are 
dealing with a symmetric one-dimensional problem we need only consider x > 
in the pre-focusing regime. This allows us to exploit the explicitly known 
Barenblatt point mass solution. A similar computation could be carried out if 
there were other explicit focusing solutions. 
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